Complex dynamics and scale invariance of one-dimensional memristive networks 
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Memristive systems, namely resistive systems with memory, are attracting considerable attention 
due to their ubiquity in several phenomena and technological applications. Here, we show that 
even the simplest one-dimensional network formed by the most common memristive elements with 
voltage threshold bears non-trivial physical properties. In particular, by taking into account the 
single element variability we find i) dynamical acceleration and slowing down of the total resistance 
in adiabatic processes, it) dependence of the final state on the history of the input signal with same 
initial conditions, Hi) existence of switching avalanches in memristive ladders, and iv) independence 
of the dynamics voltage threshold with respect to the number of memristive elements in the network 
(scale invariance). An important criterion for this scale invariance is the presence of memristive 
systems with very small threshold voltages in the ensemble. These results elucidate the role of 
memory in complex networks and are relevant to technological applications of these systems. 

PACS numbers: 



I. INTRODUCTION 

Resistors whose resistance depends on the past evolu- 
tion of the system are quite common in both basic and 
applied science, and they have been known for at least 
several decades [T]. Recently, they have attracted con- 
siderable interest in the context of memory applications 
- where they are oftentimes referred to as memristive 
systems [2] - but their range of applicability spans var- 
ious disciplines as diverse as non-traditional computing 
and biophysics [3H5]. In addition to their ubiquity, their 
theoretical description is quite simple: if a memristive 
system is subjected to a voltage V(t), its resistance can 
be written as R(x,V, t), namely it may depend on the 
voltage itself (which would simply make it a non-linear el- 
ement), and, most importantly, it depends on some state 
variable(s), x, which could be, e.g., the spin polarization 
[51 [7] or the position of atomic defects [5] , or any other 
physical property of the system that gives it memory. 

While most of the research so far has focused on the 
properties of these single elements - both by identifying 
the physical mechanisms for memory, and by devising 
innovative ways to employ them in practice - very little 
is known about their response when they are organized 
into networks, with the associated (and inevitable) ele- 
ment variability. In other words, the statistical properties 
of networks of memristive systems are largely unknown 
and, as we demonstrate below, are different than those 
of traditionally studied networks [SHU] ■ 

Addressing this issue has several immediate benefits. It 
is not at all obvious how a network whose elements have 



memory of past dynamics evolves collectively in time, or 
whether it possesses fundamental (and universal) charac- 
teristics that can be found in the natural world. In ad- 
dition, a resistance with memory is just a particular case 
of a general class of response functions that depend on 
state variables |12| - memcapacitive and meminductive 
[T2"j systems are similarly defined. Moreover, networks 
with memory have been shown to be promising candi- 
dates for the solution of complex optimization problems 
(see, e.g., Ref. [15} , and their circuit applications are 
likely to involve a combination of several elements with 
memory. Finally, it is worth mentioning that the human 
(and animal) brain is - in its most basic description - sim- 
ply a network with memory. This type of research may 
thus have consequences in neuroscience. 

With these motivations in mind, in this paper we set to 
study, both numerically and theoretically, the statistical 
properties of the simplest network of memristive systems: 
a finite ID network (Fig. [T|. In order to obtain its gen- 
eral features, we abstract as much as possible, and do not 
consider any particular physical mechanism for memory. 
However, due to physical constraints, memristive systems 
are commonly found with a threshold voltage, namely 
it takes a minimal voltage V t to change their resistance 
(see, e.g., Refs. [Tl H1H11)]) . Therefore, in order to be as 
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FIG. 1: (color online). One-dimensional network of N ran- 
domly oriented memristive systems Mi connected to a voltage 
source V. 
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close as possible to experimental realizations, we focus on 
such type of elements. Moreover, it is generally observed 
that the memristance R changes between two limiting 
values, R on and R ff- This experimental observation 
is also taken into account in our model. An extensive 
list of bipolar memristive devices with threshold can be 
found in the recent comprehensive review paper [T| . Spe- 
cific examples of memristive systems with threshold that 
can be used to verify our predictions include different 
metal/oxide/metal memristive nanodevices [8l Il4|. IT5] . 

Probability distribution functions /y t (V"), fn(R), 
/r dii (U), fn a ,,(R) satisfying the normalization condi- 
tions 



/ fi(y)dy = i (l) 

Jo 



where i = Vt, R, R on or R ff are used to describe thresh- 
old voltages, memristances and limiting values of memris- 
tances in the network, respectively. To understand the 
meaning of these distribution functions, we note that, 
e.g., fv t (V)dV is the probability to find a memristive 
device in the network with V t in the interval between 
V and V + dV. All other distribution functions are in- 
troduced in a similar way. Moreover, we suppose that 
systems with required fv t (V), fn on (R) and fn o ,JR) can 
be fabricated "on demand". (For example, one can ex- 
tract a subset of devices from an ensemble with Gaussian 
distribution to form a different distribution.) Even in an 
ensemble of devices fabricated with specific values of Vt, 
R on or R ff, an inevitable fluctuation of Vt, R on and 
R f f around their average values can also be described by 
the above mentioned distribution functions. Finally, we 
note that the memristance distribution function fn(R) is 
a time-dependent function, unlike fv t {Y), fR on (R) an d 
fRoss (R) distribution functions. The distribution func- 
tion fn(R) can be easily pre-set to the required form by 
applying appropriate pulse sequences to individual mem- 
ristive devices |17| . Below, we give several examples of 
the evolution of fn(R). 

In what follows, we will show that despite their appar- 
ent simplicity, memristive networks indeed show a quite 
complex dynamical behavior, and some general scale- 
invariant properties that can be readily verified exper- 
imentally. This paper is organized as follows. In Sec. [IT] 
we consider adiabatic processes and introduce two types 
of memristive dynamics. Sec. |III| studies the dependence 
of the final state on the history of applied voltage. Scale- 
invariance properties of memristive networks are demon- 
strated in Sec. IIVI In Sec. |V]we consider a memristive 
ladder and demonstrate the existence of avalanches in its 
switching dynamics. Finally, in Sec. |VI| we report our 
conclusions. 
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FIG. 2: (color online). Total memristance R tot in a net- 
work of 10 randomly oriented devices as a function of a 
slowly ramped voltage V. In this simulation, we have used 
fixed values of R on — 0.05i?o and R ff = 1.95i?o, and 
(3 = 2R / (V t max ■ s). The threshold voltages and initial mem- 
ristances of memristive systems have been randomly selected 
in the interval [0,F t max ] and [R on ,Rof /], respectively. 



II. ACCELERATION AND SLOWING-DOWN 
OF THE RESISTANCE SWITCHING IN 
ADIABATIC PROCESSES 

Let us first discuss the switching dynamics in a ID 
chain of memristive systems. For the sake of definiteness 
let us consider a simplified version of a generic threshold 
voltage model of memristive system [5] 

I = (2) 
x = p (Vm - 0.5 [\V M + V t \- \V M - V t \\) 

x [6 (V M ) 9 (R off -x)+6 (—Vm) (x - R on )] (3) 

where / and Vm are the current through and the volt- 
age across the system, respectively, x is the internal state 
variable playing the role of memristance, x = R, 9(-) is 
the step function, /3 is a positive constant characteriz- 
ing the rate of memristance change when \Vm\ > V tl V t 
is the threshold voltage, and R on and R jj are limiting 
values of the memristance R. In Eq. p), the role of 
^-functions is to confine the memristance change to the 
interval between R on and R ff- In our numerical sim- 
ulations, the value of x — R is monitored at each time 
step and in the situations when x < R on or x > R ff, 
it is set equal to R on or Roff, respectively. Importantly, 
the model given by Eqs. Q and ^ describes a system 
whose memristance can only change when the applied 
voltage magnitude exceeds the threshold voltage. In ad- 
dition, according to the sign of the applied voltage its 
resistance may increase or decrease, thus representing an 
asymmetric system, as indicated by the black thick line 
in the memristor symbols in Fig. [l] 

Fig. [2] presents the total memristance of a network of 
such elements as a function of a slowly ramped voltage for 
a particular network realization. The initial resistances 
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and the threshold voltages have been chosen randomly 
within uniform distributions. The polarity of the mem- 
ristive elements in the network has also been randomly 
selected. The voltage is then increased following an adi- 
abatic process: its increase rate is so slow that at any 
moment of time the network is in its equilibrium state. 
However, even in this simple adiabatic switch-on, one can 
easily notice two distinctive types of switching behavior. 
An abrupt - accelerated - switching, with well defined 
steps (see Fig. |2]), occurs when the memristance Rj of 
the j-th memristor in the network increases at the given 
voltage polarity. In fact, as soon as the voltage drop 
across such element exceeds its threshold voltage, the in- 
crease in memristance increases the voltage drop across 
it, thus accelerating the switching. Instead, for memris- 
tive systems connected with opposite polarity, a decrease 
in memristance decreases the voltage drop across the sys- 
tem, thus decelerating the switching. This effect is clearly 
seen in long switching tails. 

The dependence of the switching rate on the device 
polarity in 2D networks of ideal memristors |18| was pre- 
viously reported in Ref. [19 . This previous result, how- 
ever, does not take into account the threshold-type dy- 
namics of realistic memristive devices [T] . Consequently, 
in adiabatic experiments, instead of sharp steps and long 
tails predicted by us, this previous work suggests unre- 
alistic switching of memristors into their limiting states 
already at very small applied voltages and the absence of 
any subsequent evolution of R tot as shown in Fig. [2] 

The accelerated/decelerated switching behavior de- 
scribed above reflects also on the evolution of the resis- 
tance distribution function, fn(x), of the network. This 
is shown in Fig. [3] at several moments of time when 
the applied voltage Vq adiabatically ramps up from to 
\Qymax j n f ac ^ memristance distribution function 
dynamics (Fig. |3| indicates a faster switching of mem- 
ristive systems with initially higher mcmristances (thus 
with voltage drops exceeding their threshold voltages). 
Moreover, while the switching into the "off" state is al- 
ways complete (since it is accelerated), memristive sys- 
tems connected with opposite polarity may only partially 
switch toward the "on" state because of the decelerat- 
ing switching effect discussed above. Higher values of 
applied voltage would further promote their switching 
as well as initiate switching of memristive systems with 
higher Vt and lower initial resistances, requiring higher 
voltage drops to induce their dynamics. 

Fig. [3]presents the evolution of the resistance distribu- 
tion function in the case of normally distributed thresh- 
old voltages and limiting values of mcmristances. The 
overall behavior of fn{x) is very similar to that shown 
in Fig. [3] for the case of uniformly distributed device 
parameters. Specifically, in both cases, memristive de- 
vices switch into the limiting R a f f states faster than into 
R on ones. Moreover, the dynamics of total memristance, 
Rtotit) in the network with normally distributed device 
characteristics has exactly the same features as Rtot(t) 
depicted in Fig. [2] (acceleration and slowing-down of the 
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FIG. 3: (color online). Evolution of memristance distribution 
function /_r in a network of 10 randomly oriented memristive 
systems subjected to an adiabatically switched voltage Vb- In 
this simulation, we have used fixed values of R on = 0.5i?o and 
Raff = l.5R , /3 = 2R /(V t max -s), the threshold voltages 
and initial memristances of memristive systems have been 
randomly selected in the interval [0,V7™ aa: ] and [R on ,R j /], 
respectively. The memristance distribution function fn has 
been obtained by averaging over 10 4 random realizations of 
the network. 



switching ). 



III. DIFFERENT FINAL STATES WITH SAME 
INITIAL CONDITION 

Another consequence of memory is that the memris- 
tance at the final moment of time is determined both 
by the initial conditions and by the actual shape of the 
applied signal, so that if the initial conditions are the 
same, the final memristance could be different according 
to the way the input signal is applied. This effect is even 
more pronounced in memristive networks, where all el- 
ements are collectively coupled by current conservation. 
As an example we compare the evolution of memristances 
in a network subject to both a slowly ramped-up voltage 
(from to !0V t max ) and a step voltage of l0V t max applied 
at t = 0. The initial conditions (initial parameters of 
memristive systems in the network) are chosen to be the 
same. Fig. [5]shows results of our simulations for a partic- 
ular network realization of ten elements. One can clearly 
notice that in this particular situation the main differ- 
ence of final states is associated with the fourth mem- 
ristor from the left. While in the case of slowly ramped 
voltage R 4 (t = tf ina i) = i? 4 (i = 0), the sudden voltage 
pulse switches the fourth memristive element into the 
"off" state R ff- According to our observations, such a 
significant difference occurs in ~ 8% of network realiza- 
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FIG. 4: (color online). Evolution of memristance distribu- 
tion function fn in a network of 10 randomly oriented mem- 
ristive systems subjected to an adiabatically switched voltage 
Vo- In this simulation, we have used normally distributed 
Ron, Roff, and Vt with means 0.5i?o, 1-5 Ro and Vt and 
standard deviations 0.02i?o, 0.02J?o and 0.05V t °, respectively, 
/3 = 2Ro/ (Vt ■ s). The memristance distribution function fn 
has been obtained by averaging over 10 random realizations 
of the network. 




FIG. 5: (color online). Initial (top) and final (bottom) mem- 
ristances in a network of 10 memristive elements (denoted by 
the index i) subjected to a slow ramp-up and a step-like ap- 
plied voltage so that in both cases V(t = t final) ~ 10V t ma:c 
(the respective voltage shapes are indicated in the middle). 
In this simulation, we have used fixed values of R on — 0.05i?o 
and R f f = L95.Ro, P = 2R /(V t max • s), the threshold volt- 
ages and initial memristances have been randomly selected in 
the interval [0,V t maa: ] and [0.05i? ,l-957? ], respectively, and 
the polarity of the memristive systems in the network has 
also been randomly selected. However, both simulations have 
been performed starting with the same initial randomly se- 
lected network realization. 



tions at the selected simulation conditions. 

In order to understand the signal-dependent switching 
in the memristive networks, let us consider the simplest 
nontrivial case of a network of two memristive elements, 
Mi and M 2 , connected in series (see Fig. [T| . According 
to Eqs. (|2j, (|3j, such a network, subjected to a non- 
negative voltage V > 0, is described by the following 
equations 



dR l 
~dT 



dR 2 
dt 



V — Vti ) , when 



Ri 
R1+R2 



V > V n , 



and R\ < R a ff for sign 
or Ri > R on for sign "- 
0,in all other cases; 

and i?2 < Roff for sign " 
or R 2 > Ron for sign 
0, in all other cases, 



(4) 



,when ^V>V t2 , 



(5) 



where R\ and i?2, and Vti and V t 2 are the memristances 
and the threshold voltages of the first and second mem- 
ristive systems, respectively. The signs " ± " in the RHS 
of Eqs. Q, ([5]) take into account different possible po- 
larities of the memristive systems in the network. Here 
we suppose that the parameter j3 determining the switch- 
ing rate is positive and of the same magnitude for both 
memristive systems. Eqs. Q, (|5| determine a system of 
nonlinear differential equations, which must be supplied 



with the initial conditions for memristances R\ and R2 
at t = 0. 

Let us consider a situation when a constant voltage 
V is applied at t = 0. In this case, Eqs. Q, (|5| can 
be solved analytically. However, due to rather complex 
functional dependencies in the RHS of these equations, 
there lot of cases depending on the mutual relations 
between the initial memristances, threshold voltages, and 
applied voltage V. Below, we investigate several possible 
scenarios for signal-dependent switching. 

According to Eqs. Q, ([5|, any changes in mem- 
ristances occur when the voltage drop across at least 
one memristive system exceeds its threshold voltage. 
Let us consider then the case when the voltage drop 
V\ = R\V/(R\ + R2) across Mi exceeds Vti, but the 
voltage drop V2 = R2V/(Ri + R2) across M2 is smaller 
then Vt2- We also suppose that R on < Ri < Roff- In 
this case (see Eq. (|5|), the memristance R 2 does not 
change, R2 = const. Eq. Q can be integrated resulting 
in the following algebraic equation for Ri(t): 



Ri 



VR 2 

v-v n 



In 



(y-Vti)Ri-VtiRa 
(V-V t i)Ri(0)-V t iR 2 _ 

R 1 (0)±(3{V-V tl )t. 



(G) 



It should be noted that Eq. (|6| is valid while the following 
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inequalities are satisfied: 



Vi = 

v 2 



(Ri + Ri) 

R2 
(Ri + R2) 



< R„ff, 


(7) 


V>V n , 


(8) 


-v < v t2 . 


(9) 



If we have a "positive" polarity with respect to applied 
external voltage V (sign " + " before (3 in the RHS of Eq. 
Q), then as it follows from Eq. ([6| or directly from 
Eq. Q, the memristance Ri increases with time. In 
this case the voltage drop V\ increases, but V 2 decreases, 
because the total voltage drop across both memristors 
is constant, V\ + V 2 = V. Thus the memristance R\ 
increases with an acceleration (as V± increases). When 
Ri (t) reaches R f / (the corresponding moment of time 
t can be immediately calculated from Eq.(|6|), the mem- 
ristance R\ stops increasing (according to Eq. Q) and 
stays constant. 

In the case of a "negative" polarity (sign " — " before 
/3 in the RHS of Eq. @), it follows from Eq. (Js]) (or di- 
rectly from Eq. Q) that the memristance R\ decreases 
with time. In this case the voltage drop V 2 increases, and 
V\ decreases. This evolution, as determined by Eq.|6]) 
stops when memristance R\ reaches the lowest possible 
value Ron, or inequality ^ does not hold, depending 
on what happens first. If inequality (|9| becomes invalid 
first, while other inequalities Q, ^ are still valid, we 
encounter the situation when both memristances change 
(if, of course, R on < R 2 < R ff)- Let us now consider 
this case when both memristances are changed simulta- 
neously. This is possible only when the following inequal- 
ities are valid: 



Ron < Ri, 2 < Rofft 

V, = — V > Vtt 

1 (RX + R2Y tU 

i?9 

V 2 = 7 rV > V t2 . 

{R! + R 2 ) 



(10) 

(11) 

(12) 



We limit ourselves to the cases when the memristive 
elements, Mi and M 2 , are connected with the same po- 
larity. Then by taking the sum of Eqs. ^ and ([5]), and 
integrating it with respect to time we find for the total 
memristance 

JZi +R 2 = ±0(V - V n - V t2 )t + i?i(0) + R 2 (0). (13) 

By substituting this relation into Eq. Q we get a linear 
differential equation which can be easily integrated. As a 
result we find the following expressions which determine 
time dependence of memristances for both memristors: 



R 1 (t)=P(R 1 + R 2 



V-Vf, -v t . 



Vti 



R 2 (t) = -P(R 1 + R 2 ) 



Yr; 



V tl + V t 



-(R 1 + R 2 ),(U) 



-(R 1 +R 2 ),(15) 



where parameter P is determined by the initial conditions 
V a Ri(0) - V t iR 2 (0) 



P 



V tl + V t2 
From Eqs 



(i?i(0)+i? 2 (0)) v ~ v ti-v t3 .(16) 



( 13 1-( 16 1 it follows that the voltage drop V% 
increases with time for the "positive" polarity and when 
V t2 Ri(0) > V t iR 2 (0), while voltage drop V 2 decreases. 
Thus it leads to accelerating of the switching of memris- 
tor Mi and to decelerating of the switching of memristor 
M 2 . 

Having discussed the switching dynamics due to a con- 
stant applied voltage, let us compare the final states 
of memristive systems for two different bias application 
protocols. Specifically, assuming the same initial states 
of the memristive systems and the same final voltage 
magnitude, we consider the situations of the step and 
slowly ramped voltages. Our consideration is based on 
the following model parameters: i?i(0) = R 2 (0) = 0.2i?o, 
Ron = 0.05i? , R off = lMRo, V tl = V/3, V t2 = 1.1V/3, 
and the signs " + " in the RHS of Eqs. Q, 

For slowly ramped applied voltage, Mi is in the regime 
of accelerated switching (see also discussion after Eq. 
([9])). In addition, one can easily notice that the volt- 
age drop across M 2 never exceeds Vt 2 . Consequently, the 
final memristances are R\{t = t final) = Roff = l-95.Ro 



and R 2 (t — t 



final ) 



R 2 (t = 0) = 0.2i? . 



In the case of the step- like applied voltage V, with our 



specific model parameters the inequalities ( 10 1-( 12 1 are 
satisfied at the initial moment of time. As a result, both 
memristances are changing according to Eqs. (13l-(16l. 
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Because of the positivity of parameter P, the voltage V\ 
increases in time, but the voltage V 2 decreases (as it fol- 
lows from Eqs. (p~3])-( 15 1). At the same time, both mem- 
ristances i?i and i? 2 increase according to Eqs. Q, (|5j. 
At a certain moment of time the voltage drop V 2 across 
memristor M 2 becomes equal to the threshold voltage 
Vti, and further change of M 2 becomes impossible. Fur- 
ther evolution of memristance R\ is determined by Eq. 
(|6J) and lead to R\{t — tfi na i) — Roff — 1.95i?0) while 
the memristance R 2 does not change, R 2 (t = tfi na i) — 
0.33i?o- Clearly, the final state of M 2 in this case is differ- 
ent than that for the case of the slowly ramped voltage. 

We would like to emphasize that the different final 
states are obtained as a result of the collective dynam- 
ics of many memristive systems in the network. Such a 
behavior cannot be observed considering a single mem- 
ristive device described by Eqs. Q, (J3J) subjected to 
similar voltage patterns. The interaction among devices 
in the network is provided by the current that, at each 
moment of time, is determined by states of all devices in 
the network. 

Our observation of different final states can find useful 
applications in electronics, albeit in a modified form. For 
example, the sensitivity of the final states to the applied 
pulse profile can be used in signal recognition. We would 
also like to note that the response of single memristive 
devices with more complex internal degrees of freedom 
can be richer than that predicted by Eqs. pj. Even 



6 



with such single memory devices, one can expect the de- 
pendence of the final states on the voltage application 
protocol. Examples of such memory devices can be found 
in Refs. |201|2"T]. 



IV. SCALE INVARIANCE OF DYNAMICS 
THRESHOLD VOLTAGE 

In this Section, we ask the question of the probabil- 
ity of dynamics in the memristive network. Clearly, if 
a very low voltage V is applied to the network then the 
probability of dynamics is close to zero as the voltage 
across each memristive system is likely below its thresh- 
old voltage. If a high voltage is applied, the probability 
of dynamics is close to 1 as it is highly probable that the 
voltage across at least one memristive system exceeds its 
threshold. To quantify this effect, let us find a dynamics 
voltage threshold V* , which is the magnitude of the ap- 
plied voltage when, with a probability of 1/2, the state of 
at least a single memristive system will change. Formally, 
it can be found from 

p d (V*) = 1/2, (17) 
where pd is a probability of dynamics in the network. Al- 



ternatively, Eq. (17 1 can be rewritten as 1 — PdiV*) 



1/2, where the left-hand-side is the probability of no dy- 
namics. The latter occurs when the voltage across each 
memristive system in the network, Vi, does not exceed 
its threshold voltage, Vt.i- Consequently, the probabil- 
ity of no dynamics can be written as a product of single 
element probabilities 

p(Vi < V ttl )p(V 2 < V t , 2 )...p{V N < V t , N ) = 1/2, (18) 
where V t = RiV*/R to t and R tot = E?=i^- Notic 



ing that the probabilities in Eq. ( 18 1 can be expressed 



through the threshold voltage distribution function fy t 

as 



(19) 



p(y i <V tr ,i) = J fvMdy, 

Vi 

we finally arrive at the following equation 



fv t (yi)dyi j fv t (y2)dy 2 - J fv t {VN)dy N = 1/2 
Vi v 2 V N 



determining the dynamics threshold voltage V* 



(20) 



A. General case, narrow distribution of initial 
memristances 

Assuming that initial memristances are narrow dis- 
tributed around a certain value R, we use Vi = V* /N 



in Eq. ( 20 ) to derive 



fv t {x)dx 



V/N 



N 



= 1/2. 



(21) 



Let us then find the dynamics threshold voltage V* in the 
limit of large N. By using the normalization condition, 
and by taking the logarithm of both sides of Eq. (21 ) we 
find, for arbitrary N, the following relation 



In 



V/N 



f Vt (x)dx 



In 2 



(22) 



From Eq 
in Eq 



(22 1 we conclude that for large N the integral 
22| tends to zero. Because of the positivity of 



the distribution function fy t (x) this means that the inte- 
gration domain must tend to zero as N tends to infinity. 
Thus we see that limjv^oo V* /N = 0. This means that 
we can use the asymptotic behavior of the distribution 
function to find the dynamics threshold voltage. Let us 
then consider a general situation when 



/vt(x) = / x 7 , x^O, 



(23) 



with some positive parameters /o and 7 > 0. In this case 
the integration in Eq. ( 22 1 gives the final result 



V* = 



(1 + 7) In 2 
fo 



1/(1+7) 



N -> 00. 



(24) 



When there is a finite probability of arbitrarily small 
thresholds (7 = and f 7^ 0) , we obtain a scale invariant 
expression 



V* = 



In 2 



N -> 00. 



fvAoy 

This is the main result of Section II VI 



(25) 



B. Exponential distribution of threshold voltages 

Note that the condition Ri = R is not always nec- 
essarily required to obtain Eq. ( [25| . To see this, let 
us consider, e.g., an exponential distribution of thresh- 
old voltages, fv t {y) — exp (— y/Vo) /Vq, where Vq is a 
characteristic threshold voltage of memristive systems. 
It follows from Eq. p0| that 



_V L _V2 

e v o e v o 



1/2. 



(26) 



Since Y] Vj = V* , we find the following scale invariant 
expression for the dynamics threshold voltage: 



V* = ln(2)Vb, 



(27) 



which is in full agreement with Eq. (251. Note, however 



that Eq. ( 27 1 has been obtained without the assumption 



of narrow distribution of memristances. 
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FIG. 6: (color online). Dynamics threshold voltage as a 
function of the number of memristive systems in the network. 
This plot has been obtained assuming uniform probability 
distributions of threshold voltages and initial memristances. 
Specifically, initial memristances are selected randomly in the 
interval between Ro — AR and Ro + AR, and threshold volt- 
ages are selected in the interval between and Vm ■ The solid 
line was obtained using Eq. ( 28 1 . 



C. Uniform distribution of threshold voltages 

Let us now consider the uniform distribution of the 
threshold voltages, namely, fv{ x ) = 1/Vm, if < x < 
Vm, and fv{x) = otherwise. Here, Vm is the max- 
imum threshold voltage. For the sake of simplicity, we 
also assume a narrow distribution of initial memristances 
Ri — R. Then, from Eq. (20), it can be found that 



V* = NV M 1 - 



JV- 



-> ln(2)V 



M ■ 



(28) 



Clearly, the right hand side of Eq. ( 28 1 is again of the 
form of Eq. (25 1. Numerical calculations shown in Fig. 



[6] indicate that a random uniform distribution of initial 
memristances does not change the result predicted by Eq. 
p8l. 



D. xexp(— x) distribution of threshold voltages 

Finally, we consider an ensemble of memristive systems 
described by the threshold voltage distribution function 
fv(x) = (x/V u 2 )exp(— x/Vq). The difference with the 
previously considered cases is that now /y(0) = and 
Eq. ( |25[ ) is not applicable. Assuming a narrow distribu- 
tion of initial memristances Ri = R, Eq. ( 20 1 is written 
as 



V* 
NVa 



N 



1/2. 



(29) 



Fig. [7] shows a numerical solution of Eq. (29 1. The 
asymptotic behavior of the dynamics threshold voltage 




FIG. 7: (color online). Dynamics threshold voltage as a 
function of the number of memristive systems in the network. 
This plot has been obtained as a solution of Eq. ( 29 1 based 
on xexp(-x) distribution of threshold voltages. 



(24) with 
V \ / 2hi2N as TV -> 00. 



V* as N tends to infinity is described by Eq 
7=1, and f = l/V 2 , i.e. V 
Of course, the same result can be derived from Eq. (29 1 
as well. It is clearly seen that the dynamic threshold 
voltage V* does not saturate as in the case of previously 
considered distributions (see Eqs. (25 1, (27) and (28)). 



We relate this observation to the absence of memristive 
systems with zero V tr in the ensemble. 



V. SWITCHING AVALANCHES IN 
MEMRISTIVE LADDERS 

Finally, let us consider the possibility of avalanches in 
memristive networks. By an avalanche, we mean the situ- 
ation where a single rapid switching event induces switch- 
ings in other memristive systems. It follows from our pre- 
vious considerations (reported in Sec. [TlJ) that avalanches 
are not possible in purely ID networks. Indeed, an ac- 
celerated switching can not induce an avalanche since 
the increase in memristance of the switching element re- 
duces the voltages across all other elements in the net- 
work. During the slowing-down switching, the inducing 
switching event is not well characterized as being spread 
in time. 

Avalanches, however, are possible in networks of higher 
dimensions. A memristive ladder (a quasi-one dimen- 
sional memristive network) presented in Fig. |8ja) is an 
example of such situation. We have performed a nu- 
merical simulation of the dynamics of memristive lad- 
der consisting of 2 chains of 8 memristive systems each 
(16 memristive devices in total) subjected to a slowly 
ramped voltage. In the particular realization of memris- 
tive devices whose evolution is presented in Fig. [8^b), 
the simultaneous switching of 3-rd and 4-th memristive 
systems into the "off" state is an avalanche. 

In the memristive ladder, an avalanche is possible in 
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FIG. 8: (color online), (a) Memristive ladder consisting of N 
randomly oriented memristive systems, (b) Individual mem- 
ristances Ri in a memristive ladder of 16 randomly oriented 
devices as a function of a slowly ramped voltage. The mem- 
ristance of each device is represented by a color in the band 
in the vertical direction representing time evolution. In this 
simulation, we have used fixed values of R on = 0.5i?o and 
Raff = 1.5.Ro, P = 2R /{V t max ■ s). The threshold voltages 
and initial memristances of memristive systems have been 
randomly selected in the interval [O,^"" 131 ] and [R on ,R f f], 
respectively. 



a pair of memristive systems connected in parallel (Mj 
and Mj+i, where j is an odd number) that are in the 
regime of accelerated switching. The switching of the 
memristive system with a lower V t , say Mj, initiates the 
avalanche (the switching of Mj +1 ) if the voltage across 
these two memristive systems after the switching of Mj 
exceeds the threshold voltage of Mj+i. As the switch- 
ing in the accelerated switching regime occurs fast, effec- 



tively, the switchings of Mj and Mj +1 are observed at the 
same value of slowly ramped voltage (see Fig. [8](b) for an 
example). We expect that the size and role of switching 
avalanches increase with the increase of network dimen- 
sionality. Detailed studies of avalanches in 2D and 3D 
cases will be reported elsewhere. 



VI. CONCLUSIONS 

In conclusion, we have found that one-dimensional 
memristive networks exhibit several distinctive features 
that make them quite unlike traditionally studied dynam- 
ical networks |10[ HT] . The essential difference is related 
to intrinsic non-local interactions coupled to memory fea- 
tures. Indeed, the order of memristive systems in ID net- 
works is completely irrelevant since the current through 
the network is conserved. Two possible types of switch- 
ing dynamics - accelerating and decelerating - have been 
identified depending on the polarity of memristive ele- 
ments in the network. We have also demonstrated that 
the final state of the network depends on the protocol 
of how the external voltage is applied, and that switch- 
ing avalanches can occur in memristive ladders. Addi- 
tionally, a scale invariance of dynamics threshold voltage 
has been demonstrated. Most of our results are univer- 
sal in the sense that essentially they do not depend on 
the selected distribution functions (e.g., the acceleration 
and slowing down of switching, different final states, ex- 
istence of switching avalanches). The results reported in 
this paper are also relevant to networks of memcapacitive 
and meminductive systems [IS], while specific features of 
switching dynamics with such elements may not coin- 
cide. In general, we expect some of these predictions to 
be shared by arbitrary complex networks with memory. 
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